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Abstract. More than 600 high resolution spectra of stars with spectral type F and later were obtained in order to search for 
signatures of differential rotation in line profiles. In 147 stars the rotation law could be measured, with 28 of them found 
(f) ' to be differentially rotating. Comparison to rotation laws in stars of spectral type A reveals that differential rotation sets in 

at the convection boundary in the HR-diagram; no star that is significantly hotter than the convection boundary exhibits the 
signatures of differential rotation. Four late A-/early F-type stars close to the convection boundary and at v sini « lOOkms -1 
show extraordinarily strong absolute shear at short rotation periods around one day. It is suggested that this is due to their 
small convection zone depth and that it is connected to a narrow range in surface velocity; the four stars are very similar in T bB 
and v sin i. Detection frequencies of differential rotation a = AQ./SI > were analyzed in stars with varying temperature and 
rotation velocity. Measurable differential rotation is more frequent in late-type stars and slow rotators. The strength of absolute 
shear, An, and differential rotation a are examined as functions of the stellar effective temperature and rotation period. The 
ly-^ ' highest values of AQ are found at rotation periods between two and three days. In slower rotators, the strongest absolute shear 

at a given rotation rate Afi max is given approximately by Afi max oc /> , i.e., a max » const. In faster rotators, both a mwt and Afi max 
diminish less rapidly. A comparison with differential rotation measurements in stars of later spectral type shows that F-stars 
exhibit stronger shear than cooler stars do and the upper boundary in absolute shear Afl with temperature is consistent with the 
temperature-scaling law found in Doppler Imaging measurements. 
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1. Introduction magnetic activity. While early magnetic braking may be due 

to fossil fields amplified during contraction of the protostel- 

t3 Stellar rotation rates range from those too slow to be detected , , , .. .. . , . , . , 

CZ b lar cloud, magnetic activity at later phases requires a mecha- 

by Doppler broadening up to rates at which centrifugal forces . . . . ,. c , .. , T 

J rr b r b nism that maintains magnetic fields on longer timescales. In 

become comparable to surface gravity. Surface magnetic fields, c .. , , . , , . , t , . 

r b 1 b the Sun, a magnetic dynamo located at the interface between 

either fossil or generated by some type of magnetic dynamo, .. , , ,. .. . , . , ,. , ,. £ 

b 1 Jtf . convective envelope and radiative core is driven by radial dif - 

can couple to ionized plasma and brake a star's rotation via a t t . , . t . r„_- . , . , 

r r ferential rotation. This dynamo has been identmed as the main 

magneto-thermal wind. Magnetic braking is observed in stars t .. c . ■ ■ ■ , 

b b b source of magnetic fields, although there is growing evidence 

with deep convective envelopes where magnetic dynamo pro- . .. . . , t .. c , 

r r b 1 r that it is not the only source of magnetic held generation (e.g., 

cesses can efficiently maintain strong magnetic fields. As aeon- |„ , . — T~Z 1 Urvrvnk j .i. . - c u i 

J i,i Schrnver & Zwaan, 2000) and that magnetic helds also exist in 
sequence, field stars of a spectral type later than F are generally fuUy convective stars J.Tohns-Krull & Valentii ri996). The solar- 
slow rotators with surface velocities below 10km s . . i , i .i . .. i . 

type dynamo, however, has the potential to generate magnetic 

Magnetic braking requires the existence of a magnetic field, fiddg ^ ^ non . fully convect ive stars, as long as they show a 

which also causes the ple thora of all the effects found in stellar convective envelope, and magnetic fields have b een observed 

* Based on observations carried out at the European Southern in a variety of slowly rotating stars (e.g.,|Marcyl[IlsilQiail 

Observatory, Paranal and La Silla, PIDs 68.D-0181, 69.D-0015, 71.D- ll984ISolankill99lb. 
0127, 72.D-0159, 73.D-0139, and on observations collected at the 

Centro Astronomico Hispano Aleman (CAHA) at Calar Alto, op- Calculations of stellar rotation laws, which describe an- 

erated jointly by the Max-Planck Institut fur Astronomie and the gular velocity as a fu nction of radius and l a titude, have 

k if/inotinnu Xt Dii/iKtar 



Instituto de Astrofisica de Andalucia (CSIC) been carried out by iKitchatinov & Rudiger (1999) and 

** Marie Curie International Outgoing Fellow Kiiker & Rudiger ( 2005 ) for different equatorial angular veloc- 
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ities. iKitchatinov & Riidigerl dl999l) investigated rotation laws 
in a G2- and a K5-dwarf. They expect stronger latitudinal dif- 
ferential rotation in slower rotators, and their G2-dwarf model 
exhibits stronger d ifferential rotation than the K5 -dwarf does. 
iKiiker & Riidigerl J2005h calculate a solar-like model, as well 
as a model of an F8 main sequence star. They also come to 
the conclusion that differential rotation is stronger in the hot- 
ter model; the maximum differential rotation in the F8 star is 
roughly twice as strong as in the G2 star for the same viscosity 
parameter. The calculated dependence of horizontal shear on 
rotation rate, however, does not show a monotonic slope but has 
a maximum that occurs near 10 d in the F8 type star and around 
25 d in the solar-type star. The strength of differential rotation 
depends on the choice of the viscosity parameter, which is not 
well constrained, but the trends in temperature and rotation are 
unaffected by that choice. 

Observational confirmation of solar radial and latitudinal 
differential rotation comes from helioseismological studies that 
provide a detailed pict ure of the d ifferentially rotating outer 
convection zone (e.g.. lSchoulll998l) . Such seismological stud- 
ies are not yet available for any other star. Asteroseismological 
missions, like COROT and Kepler, may open a new window 
on stellar differential rotation, but its data q uality may provide 
only a very limited picture in the near future jGizon & Sol anki. 
2004j). In the case of the Sun, radial differential rotation man- 
ifests itself in latitudinal differential rotation that can be ob- 
served at the stellar surface, but all stars except the Sun are at 
distances where their surfaces cannot be adequately resolved. 
With the advent of large optical interferometers, more may be 
learn ed from observations of spati ally resolved stellar surfaces 
tDomiciano de Souza et all 12004). For now, we have to rely on 
indirect methods to measure the stellar rotation law. 

Photometric programs that search for stellar differential ro- 
tation assume that starspots emerge at various l atitud e s with 

different rotation rates, as observed on the Sun. Hall (1991) 

I 1 1 ' 1 — 

and Donahue et al. ( 1996) measured photometric rotation peri- 
ods and interpreted seasonal variations as the effect of differen- 
tial rotation on migrating spots. Although these techniques are 
comparable to the successful measurements of so lar differential 
rotation through sunspots (Baltha sar et all ll986). they still rely 
on a number of assumptions, e.g., the spot lifetime being longer 
than the observational seque nce, an assumption difficult to test 
in stars other than the Sun (cf. lWolter et alll2005l) . Photometric 
measurements report lower limits for differential rotation on the 
order of 10% of the rotation velocity (i.e. the equator rotating 
10% faster than the polar regions), and differential rotation is 
reported to be stronger in slower rotators. 

Doppler Imaging (DI) has been extensively used to deter- 
mine latitudinal differential rotation where the derived maps 
are constructed from time-series of high-resolution spectra. 
Differential rotation can then be extracted from comparing two 
surface maps taken with time separation of a few rot ation 
period s (e.g., bonati & Collier Cameron! 1 19971: IWolter et all 
2005). Doppler maps can also be constructed that i ncorporate 
differential rotation during the inversion algorithm dPetit et all 
120021) . Detections of diff erential rotation thro ugh DI have re- 
cently been compiled by Barnes et al . (2005), who also ana- 
lyze dependence on stellar rotation and temperature, and then 



compare them to results obtained from other techniques. Their 
results will be discussed in Sect.[6] 

The technique employed in this paper is to search for lat- 
itudinal differential rotation in the shape of stellar absorption 
line profiles. This method is applicable only to stars not dom- 
inated by spots. From a single exposure, latitudinal solar-like 
differential rotation - i.e. the equator rotating faster than polar 
regions - can be derived by measuring its unambiguous finger- 
prints in the Fourier domain. The foundatio ns of the Fo urier 
transform method (FTM) were laid bvlGravMl977tll998l) . and 
first models were done bv lBninind ( 119811) . A detailed descrip- 
tion of the fingerprints of solar-like differential rotatio n and 
the first succes sful detections are given in iReiners & Sch mitt 
(2002a, 2003a). The FTM is limited to moderately rapid ro- 
tators (see Sect.|2I>, but the big advantage of this method is 
that latitudinal differential rotation can be measured from a 
single exposure. This allows the analysis of a large sample 
of stars with a comparably small amount of telescope time. 
IReiners & Schmittl ll2003ah report on differential rotation in 
ten out of a sample of 32 stars of spectral types F0-G0 with 
projected rotation velocitie s 12kms < v sin; < 50kms _1 . 
IReiners & Schmittl d2003bl) investigated a sample of 70 rapid 
rotators with v sin/ > 45kms~ 1 and found a much lower 
fraction of differential rotators. Differential rotation has also 
been sought in A-stars t hat have no deep convective envelopes. 
IReiners & Roverl ll2004) report on three objects out of 76 in the 
range A0-F1, 60kms~ I < v sin/ < 150kms _1 , which show 
signatures of differential rotation. 

In the cited works, differential rotation is investigated in 
stars of limited spectral types and rotation velocities. In this 
paper, I aim to investigate all measurements of differential ro- 
tation from FTM, add new observations, and finally compare 
them to results from DI and theoretical predictions. Currently, 
more than 600 stars were observed during the course of this 
project, and in 147 of them the rotation law could be measured 
successfully. 

2. Differential rotation in line profiles 

Latitudinal differential rotation has a characteristic fingerprint 
in the shape of the rotational broadening that appears in each 
spectral line. Since all other line broadening mechanisms, like 
turbulence, thermal, and pressure broadening, etc., also affect 
the shape of spectral lines, the effects of differential rotation are 
very subtle. Signatures of differen tial rotation and a recipe fo r 
measuring them are presented in IReiners & Schmittl (Eo02a). 
Since the signatures are so small, and spectral line blending 
is a serious issue even in stars of moderate rotation rates at 
v sin / m 20kms~ I , single lines are not measured, but instead 
a total broadening function is constructed from many lines of 
similar intrinsic shape. This process typically involves 15 lines 
in slow rotators observed at very high resolution (Sect.[3j, and 
300 lines in rapid rotators observed at lower resolution. In order 
to derive a unique broadening function at the required preci- 
sion, high data quality is required. Detaile d information about 
demands on data quality can be found in IReiners & Schmittl 
(2003aB- 
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Interpretation of the profile's shape with rotational broad- 
ening requires that the line profiles are not affected by starspots, 
stellar winds, spectroscopic multiplicity, etc.. The Fourier 
transform method is therefore limited to unspotted single stars 
with projected rotation velocities v sin/ ^ 12kms _I , and an 
upper limit at v eq ~ 200 km s _1 is set by gravitational darken- 
ing. More specific information on the influence of starspots, of 
very rapid rotation, and some examples of detected si gnatures 
of diff erent ial rotation in line profiles can be fo und in R einersl 
d2003h . and lReiners & Schmit i lT2002aU2003albl) . 

From the derived broadening profile, the rotation law is 
determined by measuring the first two zeros of the profile's 
Fourier transform, q\ and q2- In sufficiently rapid rotators, 
those are direct indicators of rotational broadening since other 
broadening mechanisms (for example turbulence or instrumen- 
tal broadening) do not show zeros at such low frequencies. This 
is also the reason only stars with v sin/ ^ 12kms _1 can be 
studied; in slower rotators, line broadening is dominated by 
turbulence and the zeros due to rotation cannot be measured. 
In stars spinning fast enough to be analyzed with FTM, ap- 
proximation of net broadening by convolutions is also justi- 
fied. The important point in choosing the Fourier domain for 
profile analysis is that convolutions become multiplications in 
Fourier space. Thus, the fingerprints of rotational broadening 
are directly visible in the observed broadening profile's Fourier 
transform, and the spectra do not have to be corrected for in- 
strumental or for any other line broadening, as long as the tar- 
gets' rotation dominates the important frequency range. 

In the following, the stellar rotation law will be approxi- 
mated in analogy to the solar case. Differential rotation is ex- 
pressed in terms of the variable a, with Q the angular velocity 
and / the latitude. The rotation law is approximated by 

Q(0 = £l Et]ualor ■ asm 2 (l), (1) 

^Equator ~ Qpole AO 

a = = . (2) 

no 

1 ^Equator 1 a 

Solar-like differential rotation is characterized by a > 
(a Q m 0.2). Sometimes AO is called absolute differential ro- 
tation. To avoid confusion with relative differential rotation 
a = AO/O, I will refer to AO as absolute shear and to or as 
differential rotation. A O is giv en in units of radd 1 . As shown 
in R einers & SchmitJ ll2003al) . the parameter a/ Vsin/, with / 
the inclination of the stellar rotation axis, can be directly ob- 
tained from the ratio of the Fourier transform's first two zeros, 
qilqx- Thus, determination of the stellar rotation law in Eq.^ 
from a single stellar spectrum is a straightforward exercise. 

2.1. Variables to describe differential rotation 

The role of differential rotation especially for magnetic field 
generation in rapid rotators is not understood well. In the case 
of the Sun, we know that its magnetic cycle is driven, at least 
in part, by radial differential rotation, which itself is reflected 
in latitudinal differential rotation. How (and if) magnetic dy- 
namos depend on the strength of differential rotation in terms of 
a or on absolute shear AO has not yet been empirically tested. 
Furthermore, different observing techniques measure different 



quantities, and authors express their results on stellar rotation 
laws in different variables. A variable frequently used is the 
lap time, i.e. the time it takes the equator to lap the polar re- 
gions (or vice versa), which is essentially the reciprocal of the 
shear AO. Differential rotation a = AO/O is the shear divided 
by the angular velocity. Expressing rotation velocity in terms 
of rotation period P essentially means the reciprocal of angular 
velocity O, but using AP instead of AO introduces several prob- 
lems, since one has to consider whether P denotes equatorial or 
polar rotation period, the latter being larger with solar-like dif- 
ferential rotation (cp. R einers & Schmittll2003bl) . In this paper, 
I will express the rotation law in terms of a and AO, in order 
to search for correlations with rotation velocity. The quantity 
measured by the FTM is a with a fixed observational threshold 
of ffmin ~ 0.05 thereby limiting the detection of deviations to 
rigid rotation. The uncertainty in a measured by FTM is ap- 
proximately 6a - 0.1. The perhaps more intuitive parameter 
for the physical consequences of differential rotation, however, 
is the shear AO = aO, where O has to be obtained from v sin i 
and the radius. When analyzing the rotation law in terms of 
AO, one has to keep in mind that it can never exceed angular 
velocity itself, if polar and equatorial regions are not allowed 
to rotate in opposite directions. In other words, differential ro- 
tation cannot be larger than 100%. 

3. Observational data 

Data for this study have been compiled from observations car- 
ried out at different telescopes. Observations of field stars with 
projected rotation velocities higher than v sin/ = 45kms _1 
have been carried out with FEROS on the 1.5m telescope 
(R = 48 000), ESO, La Silla, or with FOCES (R = 40000) 
on the 2.2m at CAHA. Slower rotating field stars were ob- 
served at higher resolution at the CES with the 3.6m telescope, 
ESO, La Silla (R = 220000). Additionally, observations in 
open clusters were carried out with the multi-object facility 
FLAMES feeding VLT's optical high-resolution spectrograph 
UVES at a resolution of R — 47 000. Details of the spectra and 
instruments used are given in Ta bleU] Parts of the data were 
published in Rein ers & SchmitJ d2003albh and iReiners et alJ 
(2005). For more information about observations and data re- 
duction, the reader is referred to these papers. Ten of the cluster 
targets (FLAMES/UVES), as well as 34 FOCES, targets were 
not reported on in former publications 1 . 

More than 600 stars were observed for this project dur- 
ing the last four years. For this analysis, I selected the 147 
stars for which the ratio of the Fourier transform's zeros q2/qi 
- the tracer of the rotation law - is measured with a preci- 
sion better than 0.1, i.e. better than 6% in case of the typ- 
ical value of qi/qi = 1.76. These 147 stars exhibit broad- 
ening functions that are (i) symmetric (to avoid contamina- 
tion by starspots), and (ii) reveal rotation velocities between 

1 It should be noted that the FLAMES/UVES targets observed in 
open cluster fields are not necessarily cluster members, as was shown 
for example for CI* IC 4665 V 102 in lReiners etafl <2005h . A detailed 
investigation of cluster membership goes beyond the scope of this 
analysis; rotation velocity and spectral type make them sufficiently 
comparable to the (probable) field stars in this context. 
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Table 1. Data used for this analysis. Spectra were taken with different high-resolution spectrographs. Detection of differential 
rotation requires higher resolution in slower rotators. The minimum rotation velocities required for determining the rotation law 
are given in column 4. From the 147 spectra from which differential rotation could be determined, 28 stars show signatures of 
differential rotation 



Instrument 


Resolving power i/A/t 


Wavelength Range [A] 


min. v sin i [km/s] 


total # of stars 


# of differential rotators 


CES 


220000 


40 


12 


32 


14 


FEROS 


48 000 


6000 


45 


23 


3 


FOCES 


40000 


6000 


45 


81 


7 


FLAMES/UVES 


47 000 


2000 


45 


11 


4 


total 








147 


28 



v sin/ = 12kms _I and 150kms~', the latter being an arbitrary 
threshold in order to minimize the amou nt of gravity-darkened 
stars with v eq » 200 km s 1 llReinersl 120031) Depending on 
the spectral resolution, the minimum rotation velocity is higher 
than v sin / = 12kms _1 (CES) or 45 km s 1 (FEROS, FOCES, 
UVES, see also TableQ). 

More than 450 stars have not been used in this analysis. 
Their broadening functions show a whole variety of broadening 
profiles. Many are slow rotators, although some were repo rted 
as rapi d rotators in earlier catalogues (cf. IReiners & Schmitl 
2003b). A large number of spectroscopic binaries (or even 
triples) were found, and another large part of the sample shows 
spectra that are apparently distorted by starspots or other mech- 
anisms that cause the spectra to appear asymmetric. For ana- 
lyzing the rotation law from the broadening profile alone, these 
spectra are useless and are not considered in this work. They 
may be promising targets for DI techniques or other science; a 
cat alogue of broaden i ng pro files is in preparation. 

Rein ers & Roven J2004I) measured the rotation law in 78 
stars of spectral type A including a few early F-type stars. 
This sample will only be incorporated in the HR-diagram in 
Sect. 16.11 As will be shown there, rotation laws in A-stars are 
fundamentally different from those in F-stars. For this reason, 
I will not incorporate the A-stars in the sample analyzed for 
rotation- and temperature-dependencies of differential rotation. 
One of the stars observed during this project has recently been 
studied with DI, with comparison of the results given in the 
Appendix in Sect.lAl 

4. Accuracy of rotation law measurements 

Projected rotation velocities v sin i were derived from the 
first zero q\ in the Fo urier transform as explained in 
IReiners & Sch mitt (2003b). As mentioned there, the precision 
of this measurement is usually < 1 kms -1 , but simulations re- 
vealed a systematic uncertainty of ~ 5% in v sin i. Thus, I chose 
the maximum of the uncertainty in the intrinsic measurement 
and the 5% limit as the error on v sin i. 

The typical uncertainty for the determined values of rela- 
tive differenti al rotation in terms of a — AQ/Q is 6a « 0.1 
( IReiners & Schmitl l2003al) . In slow rotators, this uncertainty 
is dominated by the noise level due to comparably sparse sam- 
pling of the profile even at very high resolution. In fast rotators, 
the profile (i.e., the zeros of the profile's Fourier transform) can 



be measured with very high precision. Here the uncertainty in 
a stems from poor knowledge of the limb darkening parame- 
ter e. For the case of a linear limb darkening law, the value of 
the measured ratio of the Fourier transform's first two zeros, 
Q\ an d q2, is in the range 1.72 < q2/q\ < F85 foravins et all 
1990). With the very conservative estimate that e is essentially 
unknown (0.0 < e < 1.0), every star with 1.72 < ^2/^1 < 1-85 
in this analysis was considered a rigid rotator, and qilqi < 1 .72 
was interpreted as solar-like differential rotation with the equa- 
tor rotating faster than the polar regions. From the measured 
value of qilqx, the parameter a/ V sin / was determined as ex- 
plained above. 

Stars can exhibit ratios of qi/qi > 1-85 as well, but this 
value also does not agree with rigid rotation of a homogenous 
stellar surface. In the investigated sample of 147 stars, eight 
(5%) exhibit a ratio of qilq\ > 1-85. In contrast to solar-like 
differential rotation with the equator rotating more rapidly than 
polar regions, this case may be due to anti solar-like differential 
rotation with polar regions rotating faster than the equator. On 
the other hand, it can also be caused by a co ol polar cap, w hich 
is exp ected to occur in rapidly rotating stars JSchriiver & TitleL 
2001). The lower flux emerging from a cool polar cap makes 
the line center shallow and has the same signature as anti solar- 
like differential rotation. While anti solar-like differential rota- 
tion cannot be distinguished from a cool polar cap by investi- 
gating the rotation profile, the existence of a cool polar cap in 
rapidly rotating F-stars seems much more plausible than anti 
solar-like differential rotation. Differentiation between the two 
cases can only be achieved by measuring differential rotation 
independent of the line shape. For the scope of this paper, how- 
ever, I will interpret stars with q2/qi > 1.85 as rigid rotators 
with a cool polar cap. 

5. Rotation velocity and spectral type 

The observations carried out with FEROS, FOCES, UVES, and 
the CES are very homogeneous in terms of data quality and the 
investigated wavelength region, and only measurements that 
fulfill the criterion 6a < 0.1 are considered. The targets have 
similar spectral types, but they do not form a statistically unbi- 
ased sample. The quality requirements discussed above make 
it difficult to analyze stars in a well-defined sample. This sam- 
ple certainly is severely biased by observational and systematic 
effects. The most important bias is probably due to rotational 
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Table 2. Amount of differential rotators, total number of stars, 
and percentage of differential rotators for each subsample of 
Fig-IU 



Subsample 


# of diff. rotators 


total # of stars 


percentage 


I 


2 


57 


4% 


11 


4 


24 


17% 


111 


5 


21 


24% 


IV 


17 


45 


38% 
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Fig. 1. Projected rotation velocity plotted vs. B - V for the 
whole sample of 147 stars. Differential rotators are plotted as 
filled circles. Contents of subsamples I-IV are given in Table|2] 

braking, by connecting spectral type with rotational velocity 
in the presumably old field stars. Later spectral types suffer 
from stronger magnetic braking and are expected to be gen- 
erally slower than earlier spectral types. 

Projected rotation velocities v sin; of the sample stars are 
plotted versus color B - V in Fig.^ while differentially rotat- 
ing stars (i.e. stars with a > 0) are indicated by filled circles 
and will be discussed in detail in Sect.[6] A clear dependence 
of v sin i on color is apparent, as expected. In Fig.^ the sample 
is divided arbitrarily into four subsamples with projected rota- 
tion velocities that are higher (lower) than v sin / = 70kms~ 1 , 
and color redder (bluer) than B-V - 0.4. The total numbers of 
stars in each subsample and the numbers of differential rotators 
are given in Table[2] Most targets occupy regions I and IV. The 
scarcity in region II is due to the early F-type stars not being 
subject to strong magnetic braking; since the sample mainly 
consists of field stars, most stars later than B - V - 0.4 have 
been decelerated into region IV and do not appear in region 
III. Thus, slower rotating stars generally have a later spectral 
type in the sample. This implies that the effects of temperature 
and rotation velocity on the fraction of differentially rotating 
stars are degenerate in this sample; slow rotation implies late 
spectral type. It is therefore not possible to uniquely distinguish 
between the effects of rotation and spectral type in the investi- 
gated sample. This degeneracy has to be kept in mind when 
trying to interpret rotation- and temperature-dependencies of 
differential rotation in the following chapters. 



6. Results 

As a first result of the profile analysis, I present v sin i and 
the measurement of the rotation law for all 147 objects in 
AppendixE Fields stars (FEROS, FOCES, and CES observa- 
tions) are compiled in Table lB.ll FLAMES/UVES targets in 
TablelO 

In the sample of 147 stars for which the rotation law was 
measured, 28 (19%) exhibit signatures of solar-like differen- 
tial rotation (q 2 /qi < 1.72). CES-, FEROS-, and FOCES- 
samples contain field stars of spectral type later than F0, with 
the majority of them brighter than V — 6 mag. Stars in fields 
of open clusters contained in the UVES-sample probably are 
mostly younger than the field stars, and all o f them should have 
reached the main sequence JStahler & Pallai 12004). 

6.1. Differential rotation in the HR-diagram 

All currently available measurements of the stellar rotation 
law from rotation profile shape (i.e, 147 s t ars f rom this sam- 
ple and 78 stars from iReiners & Roverl 120041) are plotted 
in an HR-diagram in Fig.|2] For the field and A-stars, effec- 
tive temperature T e ff and bolometric ma gnitude Mh n i were de- 
rived from uvbyB photometry taken from lHauck & Mermilliod l 
(Il998[) using the program UVBYB ETA published by iMoonl 
dl985h . For 7/ eff a n ew calibration bvlNapiwotzki et all dl993l) 
based on the grids of Moon & Dworetskv ( 1985) was used, and 
the statistical error of the temperature determination is about 
Ar e ff = 150 K. For three of the field stars, no uvbyB photome- 
try is available; all three are rigid rotators, so no value of T e s 
was calculated to avoid inconsistencies. For the cluster stars, 
no uvbyB data is available, so radius and temperature are esti- 
mated from B-V color using zero-ag e main sequence (ZAMS) 
polynomials taken from Gray] 1119761) . i.e., they are assumed to 
be young. For them, Mboi is calculated from 

M bol = 42.36 - 10 log r eff - 5 log/?//? , (3) 

with R e the solar radius. To get absolute V-magnitudes My 
from bolometr ic magn itudes Mboi, the bolometric correction 
calibrated from lHavesN 19781) was applied. 

In Fig. |3 stars are plotted as open or filled circles that indi- 
cate rigid (a = 0) or differential rotation {a > 0), respectively. 
Circle sizes display projected rotation velocities v sin i as ex- 
plained in the figure. The 147 stars given in TableslB~T1an dlB.2l 
are pl otted as black symbols, while stars from Reiners & Rover 
(2004) are plotted as grey symbols. Evolutionary main se- 
quence trac ks for M = 1.2, 1 .4, 1.6,2.0, and 3.0 M and the 
ZAMS from lSiess et alJ d2000l) are also shown. Near the inter- 
face of A-stars and F -stars, the "granulation boundary" from 
iGrav & Nagel dl989h is indicated with dashed lines. This is 
the region where line bisectors measured in slow rotators show 
a reversal. For dwarfs and subgiants, the "granulation bound- 
ary" coincides with theoretical calculations of the "convection 
boundary" and thus can be identifie d as the region where deep 
envelope convection disappears Ccf. lGrav & Nagell ll989). The 
exact region of this boundary is not defined well and may de- 
pend on factors other than temperature and luminosity. 
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Fig. 2. HR-diagram of all currently available measurements of latitudinal differential rotation from FTM. Stars analyzed in this 
work are plotted in black, and stars from the A-star sample in R einers & R over (2004) in grey. Differential rotators are indicated 
by filled circles. Symbol size represents ranges of projected rotation veloci ty v sin as explai ned in the figure. A typical error bar 
is given in the upper right corner. Evolutionary tr acks and the ZAMS fr om lSiess et aD fcOOOh are overplotted. Dashed lines mark 
the area of the granulation boundary according to lGrav & NageJ dl989l) . while no deep convection zones are expected on the hot 
side of the boundary. Downward arrows indicate the four stars from Table[3](see text). 



The stars shown in Fig. [2] cover a wide range in tempera- 
ture on both sides of the convection boundary. Because field 
stars are investigated, many stars have evolved away from the 
ZAMS. Rotation velocity as indicated by symbol size follows 
the well-known behavior of magnetic braking; late type stars 
are generally slowed down and exhibit slower rotation veloci- 
ties. The striking fact visible in Fig.|2]is that none of the differ- 
ential rotators detected with FTM lies significantly on the hot 
side of the convection boundary - latitudinal differential rota- 
tion has only been detected in stars believed to possess a deep 
convective envelope. 

The largest group of differential rotators can be found near 
the ZAMS at all temperatures that are cooler than the convec- 
tion boundary. The convection boundary itself is also populated 
by differential rotators. A few others can be found far away 
from the ZAMS (at masses <: 1 .6 M ) to the right end of the 
main sequence tracks. There is a hint that differential rotators 
are lacking in the region between these few rotators and those 
near the ZAMS and again between them and those near the con- 
vection boundary, although this region is occupied by rigidly 



rotating stars (i.e., stars with differential rotation weaker than 
the observational threshold). However, it is not clear from the 
available sample whether temperature (and evolutionary stage) 
or rotation velocity may be the important parameter in deter- 
mining differential rotation in F-type stars (see Sect.[5}. 

6.2. Differentially rotating A-stars at the convection 
boundary 

Among all the differentially rotating stars shown as filled cir- 
cles in Fig. |3 the four stars with the strongest shear A£2 are 
found very close to the convection boundary. This group of 
rapidly rotating late A-/early F-type stars is listed in Table|3f ■ 



2 It should be noted that the only star in this group, that lies on 
the cool side of the convection boundary, is CI* IC4665 V 102. 
Temperature and luminosity were estimated from B - V color using 
ZAMS polynomials. CI* IC 46 65 V 102, h owever, is not a member 
of the cluster IC4665 l Re iners et alll2005l) . hence its position in the 
HR-diagram is most uncertain. It may have a much higher absolute 
luminosity that shifts it even closer towards the convection boundary. 
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Table 3. Four stars with extreme values of differential rotation 
A£X It is speculated that these objects represent a special class 
of targets (marked with downward arrows in Fig. [5) in a narrow 
region of rotation velocity and effective temperatures near the 
convection boundary 



Star 


Sp 


v sini 


AQ sin i 




My 






[kms 1 ] 


[radd- 1 ] 


[K] 




HD 6869 


A9 


95 ± 5 1 


2.3 ± 1.5 1 


7390 


3.01 


HD 60555 


A6 


109 ± 5 1 


1.3 + 0.7' 


7145 


1.49 


HD 109238 


F0 


103 ± 4 1 


1.3 ±0.6' 


7184 


1.31 


IC4665 V102 


A9 


105 ± 12 2 


3.6±0.8 2 


7136 


2.93 



[ Reiners & RoyeJ <2004 
iReiners et al]l2005h 



All four exhibit very similar effective temperatures around 
T e g = 7200 K, putting them into the region where convec- 
tive envelopes are extremely shallow. The four stars show re- 
markably similar projected rotation velocities all within 10% of 
v sin / = lOOkms -1 , and all four have rotation periods shorter 
than P = 2d. Two of them, HD6869 and CI* IC4665 V 102, 
exhibit a shear that is as strong as AD. > 2radd~' at rotation 
periods that are shorter than one day. This contrasts the 34 other 
stars that show comparably small rotation periods but have dif- 
ferent temperatures or surface velocities. None of the 34 other 
rapid rotators exhibits a shear AQ sin i in excess of 0.7 rad d" 1 , 
i.e. a factor of three weaker than the two mentioned above, and 
only three of the 34 show AO > at all. 

The fact that the four strongest differential rotators are 
found at hot temperatures and high rotation velocities contra- 
dicts the general trend that differential rotation is more com- 
mon in slowly rotating cool stars, which will be discussed in 
the following sections. All four are very close to the convection 
boundary or even on its blue side, meaning extremely shallow 
convective envelopes. This leads to the suspicion that the mech- 
anism responsible for the strong shear is different from the one 
that drives the shear in stars with deeper convection zones. This 
is supported by the observation that all four stars also exhibit 
very similar (and comparably large) surface velocities of about 
100 km s _1 . It is thus suggested that the strong surface shear 
in the four stars in Table|3]is not comparable to the rest of the 
sample, but is facilitated by the high surface velocity in a par- 
ticularly shallow convective envelope. Such a mechanism could 
be supported by eigenmodes that are comparable to pulsational 
instabilities, for example in S Scu stars, but further investiga- 
tion is beyond the scope of this paper. 

6.3. The fractional amount of differential rotators 

In this section, the hypothesis that the fraction of differentially 
rotating stars with a > is independent of rotation velocity 
v sin i and color B — V will be tested. The actual strength of 
differential rotation is not taken into account, but will be inves- 
tigated in the following sections. 



In Fig.Q] the whole sample of 147 stars has been divided 
into the four arbitrary regions mentioned in Sect.[5] (dividing 
at v sin; = 70kms _1 and B — V — 0.4). The number of stars 
and differential rotators (a > 0, for this sample this means stars 
with a above the observational threshold of ff m i n as 0.05), and 
the percentage of differential rotators is given in Tabled The 
percentage of differentially rotating stars among slow rotators 
in regions II and IV is higher than it is among rapid rotators in 
regions I and III, respectively; a trend towards a higher fraction 
of differential rotators at slower rotation velocity is visible in 
both color regimes. The same is true for the percentage of dif- 
ferentially rotating stars among late-type stars in regions III and 
IV. A comparison to earlier type stars among regions I and II, 
yields a higher fraction of differentially rotating stars towards 
later spectral type. I tested the hypothesis that subsamples I-IV 
are drawn from the same distribution with a total mean of 19 % 
differential rotators. Samples II and III are consistent with this 
hypothesis (5 and 4 expected differential rotators, respectively). 
For samples I and IV, the hypothesis can be rejected at a 99 % 
level (99.3 %, 1 1 expected, and 99.6 %, 8 expected for samples 
I and IV, respectively). Thus, the fraction of differential rotators 
with a larger than the observational threshold of a m \„ w 0.05 - 
which does not depend on v sin i or color - is not constant. It is 
larger in slower rotators and stars of a later spectral type. 

The distribution of differential rotators in v sin i and B - V 
is investigated further in Fig.[5] The upper panel shows the to- 
tal number of stars divided into five bins in B - V (left), and 
seven bins in v sin i (right). The lower panel of Fig.[3]displays 
the fraction of differential rotators (a ^ 0.05) in the respective 
color/rotation bins with 2cr-errors. For example, 48% of the 25 
stars with projected rotation velocities v sin / between 10 and 
30km s show signatures of differential rotation. Although 
some bins are sparsely populated and have large errors, the 
change from color B— V = 0.4 to B - V — 0.5 and the transition 
from slow rotators to stars with v sin / > 30 km s are signif- 
icant. The trends indicated in Table|2] stand out in the lower 
panel of Fig.|5] In this sample, profiles with a ^ 0.05 are more 
frequent in slow rotators, which implies that they are more fre- 
quent in stars of later spectral type (due to the sample bias, cf. 
Sect.|5j. 

6.4. Differential rotation and rotation period 

In this and the following sections, I will investigate the strength 
of differential rotation and especially the maximum strength of 
differential rotation at different rotation rates. I discuss differ- 
ential rotation in terms of a = AQ/Q in Sect. 16.4. Tl and ana- 
lyze the shear AD in Sect. l6.4.2l Although both quantities es- 
sentially have the same meaning, it is instructive to look at both 
of them separately. Rotation speed will be expressed as a func- 
tion of equatorial rotation period instead of rotation velocity. 
The rotation period itself is not measured for the majority of 
stars, so I will use P/ sin i instead, as calculated from measured 
v sin / and the radius according to Eq.[3] 

In the left and right panels of Fig.|4] values of differen- 
tial rotation a and absolute shear AO are plotted against pro- 
jected rotation period P/ sin /, respectively. From the 147 stars 
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Fig. 3. Upper panel: Histograms of the sample in B - V (left column) and v sin / (right column). Lower panel: Fraction of stars 
with strong relative differential rotation in the bins plotted in the upper panel. 2<x-errors are overplotted. 




If- 
t I 
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Fig. 4. Differential rotation vs. projected rotation period P/ sin i. Left: Differential rotation a = AO/O; right: Absolute shear AO 
(note the logarithmic scale on the ordinate). Field F-stars are in dicated by filled c ircles, open circles display stars from cluster 
observations. Crosses m ark the three differential rotators from Reiners & Rover (2004). Short dashed lines are models from 
Kiiker & J2005I) for an F8-star with a viscosity parameter of c v = 0.15. The four stars with extraordinarily strong shear 

AO were discussed in Sect.[3] a different mechanism is suggested to cause their strong shear. For the remaining stars, the upper 
envelopes are indicated qualitatively by the long dashed lines. In the right panel, dotted lines show regions of constant differential 
rotation a — 0.1 (approximate sensitivity of FTM) and a — 1.0 (maximum if no counter-rotation of the polar regions is allowed 
for). 



in Tabled only the 28 stars with signatures of differential ro- 
tation are shown, while the other 119 objects populate the 
a = AO/O = region at 0.5 d < P/sini < ll.Od. Field 
stars from this sample with available uvbyB measurements 
are plotted as filled circles, while the four cluster targets for 
which the ZAMS-age has been assumed are plotted as open 
circl es. In both panels , the three differentially rotating stars 
from iReiners & Rover (2004), which have been discussed in 
Sect. 16.21 are marked as crosses. 

In the left and right panels of Fig.|4] the long dashed 
lines qualitatively indicate the upper envelope of a I Vsin i 
and AO Vsin i, respectively. No fit was attempted, so the lines 



should only guide the eye to clarify what will be discussed in 
the next sections. 



6.4.1. Differential rotation a 

The advantage of using a is that it is measured directly and 
its detection does not depend on rotation period, hence radius, 
while measuring AO does. Differential rotation a is smaller 
than 0.45 for all rotation periods. While a minimum threshold 
of (a/ Vsin/) m in ~ 0.05 applies, the observational technique 
has no limitations towards high values of a. Thus, the highest 
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0.45 is 



detected value of differential rotation (a/ Vsin/) m 
not due to limitations of the FTM. 

At rotation periods between two and ten days, the targets 
populate the whole region < a/ Vsin/ < 0.45, while the 
slower targets could not be analyzed due to the limitations of 
the FTM (cp. Sect[2}. Among the rapid rotators with projected 
rotation periods that are shorter than two days, the upper en- 
velope shows a clear decline among the F-stars (filled circles 
in the left panel of Fig.|4j. Except for the group of A-stars 
discussed in Sect. 16.21 (listed in Table[3j, no star with a pro- 
jected rotation period less than two days shows a > 0.2, and no 
star with Pj sin / < 1 d shows a > 0.1. Neglecting those four 
stars, the maximum value in differential rotation, (a/ Vsin/) max , 
grows from virtually zero at P — 0.5 d to (a/ Vsin /) max w 0.45 
in stars slower than P — 2 d. In slower rotators, (a/ V sin /) max 
remains approximately constant. 

6.4.2. Absolute shear AQ. 

Absolute shear (AO Vsin /) is shown in the right panel of Fig.|4] 
Since observed values of both a and O depend on inclination i, 
the observed absolute shear is 



AQ b s = ffobs^obs = a I Vsin i O sin i = AQ ' 



(4) 



In this and the following sections, I omit the factor Vsin / for 
readability. Note that in the case of small inclination angles, 
the value of AO, can be larger than AO, Vsin/, while a can 
be smaller than a/ Vsin /. In the right panel of Fig.0] the two 
dotted lines indicate the slopes of constant differential rotation 
a - AO/O. The upper line is for a — 1.0, i.e., the maximum 
differential rotation possible regardless of the observation tech- 
nique used, and the lower line shows a = 0.05, the approximate 
minimum threshold for the FTM, as explained above. 

As expected from Sect. 16.4. 11 the F-stars form a relatively 
smooth upper envelope in the maximum absolute shear AO max , 
as observed at different rotation rates. The slowest rotators ex- 
hibit low values of AQ max w 0.2radd _I at P » 10 d. AO max 
grows towards a faster rotation rate with a maximum between 
two and three days before it diminishes slightly with more 
rapid rotation. The strongest differential rotation occurs at ro- 
tation periods P between two and three days at a magnitude of 
AQ. * lradd -1 (i.e., lapping times on the order of 10 d). At 
higher rotation rates in the range 0.5 d < P < 2 d, the max- 
imum shear has a slope of roughly AQ oc P +0A . This slope, 
however, is not constrained well due to the large uncertainties 
and sparse sampling. The data are also consistent with a plateau 
at AO « 0.7radd _1 for 0.5 < P < 3 d. At slower rotation, right 
from the maximum, the slop e is approximately AO oc P l . 

Kiiker & Riidiger (2005) recently calculated AO in an F8 
star for different rotation periods. Their results for AO(P) 
(with a viscosity coefficient c v = 0.15) are displayed qualita- 
tively in both panels of Fig. El as a short-dashed line (Fig. 6 in 
Kiiker & RudigeH 1200.4 . One of their results is that in their 
model AO does not follow a single scaling relation for all pe- 
riods (as was approxim ately the case in the calculations by 
Kitchatinov & Rudigerlll999i) . but that a maximum shear arises 
at a rotation period of about P = 7 d in the case of the mod- 
eled F8-star. Comparison of their calculations (right panel of 



Fig-EJ to the upper envelope suggested in this work (long- 
dashed line) still shows a large quantitative discrepancy. The 
qualitative slopes of both curves, however, are in reasonable 
agreement with ea ch other. The theoretical curve was calcu- 
lated for an F8-star. Kiiker & Riidiger (2005) also show AO(f) 
for a solar-type star, where AO(P) is essentially moved towards 
higher rotation periods and lesser shear; i.e. the short-dashed 
curve in the right panel of Fig.|4]moves to the lower right for 
later spectr al types. Although earlier spectral types are not cal- 
culated by Kiiker & Riidiger (2005), it can be expected that 
AO(f) will shift towards higher shear and shorter rotation pe- 
riods in stars of an even earlier spectral type. Since most stars 
investigated in this sample are earlier than F8, this would sug- 
gest qualitative consistency between theoretical curves and the 
slope of AO max shown here. 

6.4.3. Comparing different techniques 

Differential rotation measurements are now available from 
a variety of observational techniques (see Sect.Q, compari- 
son of results from techniques becomes possible. However, 
such a comparison has to be carried out with great care. 
Photometrically measured periods, for example, are only sen- 
sitive to latitudes covered by spots; they reflect only parts of 
the rotation law and are always lower limits. Furthermore, tem- 
perature has been shown to be the dominating factor for the 
strength of differential rotation JKitc hatinov & Riidig erl [l 999 ; 
iKiiker & Riidiger! 120051 iBarnes et all l2005h . which has to be 
taken into account when analyzing the rotation dependence of 
differential rotation. 

In the past, analyses of relations between rotation and dif- 
ferential rotation have generally assumed a monotonic scal- 
ing relation between period P (or angular velocity O = 
2jt/P) a nd AO. Such a relation was expected from calcula- 
tions by [K itchatinov & Riidiger ( 1999|. As mentioned above, 
Kiiker & Riidiger (2005) recently presented new calculations 
showing that the AO vs. O-relation may have a temperature de- 
pendent maximum. 

Searching for dependence on angular velocitv. lBarnes etafl 
d2005h recently compiled data from differential rotation 
measurements from DI, photometric monitoring, and FTM. 
Fitting a single power law to the compiled data, they de- 
rive AO oc O +015±01 °, which is compared to the ca se of a 
G2 dwarf calculated in iKitchatinov & Riidigerl dl999l) . From 
the latter, they cite the theoretical G2 dwarf relation as 
AO oc O +0 15 , and claim agreement to t heir fit. Although 
the work of Kitchatinov & Riidiger ( 1999) has been super- 
seded by Kiiker & Riidiger (2005), it should be mentioned that 
IKitchatinov & Riidigerl i 19991) report AO oc O" 15 , implying 
stronger shear for slower rotation instead of weaker (note that 
Kiiker & Riidiger, 2005, also report a negative exponent for 
their solar-like star model at periods less than 20 d). In fact, 
the large scatter in the compilati on of all measu rements from 
different techniques (Fig. 3 in IBarnes et all 120051) and the se- 
vere bias due to systematic uncertainties (like the lack of small 
values AO in the sample measured with FTM) leads to any 
conclusion about the period dependence from such a heteroge- 
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Fig. 5. Differential rotation AQ with effective temperature T e e 
for sample stars (symbols as in Fig.0}. Big circles indicate 
large rotation velocity. Stars analyzed with DI are indicated as 
squares, and the fit to DI targets from iBarnes et alJ J2005I) is 
overplotted. F-stars with the strongest differential rotation are 
consistent with the fit. 

neous sample very uncertain. The new data in this paper does 
not significantly improve this situat ion and no analysis that im- 
proves upon the one performed by B arnes et alJ J2005h can be 
expected. 

Although the constantly growing amount of differential ro- 
tation measurements provides a relatively large sample, the re- 
sults from analyzing all measurements as one sample do not yet 
provide convincing evidence for a unique rotation dependence 
of differential rotation over the whole range of rotation periods. 

6.5. Upper bound of differential rotation and effective 
temperature 

Effective temperature is the second parameter after rotation 
to govern stellar differential rotation, since convection zone 
dept h, as well as convec tion velocity, are very sensitive to T ^. 
iKiike r & Rudiger (2005| andlKitchatinov & Rudigerldl999l) re- 
port stronger differential rotation with higher effective tempera- 
tures when comparing a solar-type star to stars of spectral types 
F8 and K5. The differential rotators in the sample of F-stars in- 
vestigated here span a range in effective temperature between 
6000 K and 7150K. As a result, analysis of temperature effects 
on differential rotation is limited by the small range of targets 
in r e ff, and is biased by the large range in rotation periods, as 
discussed above. Thus, I limit the analysis of temperature de- 
pendence to a comparison to differential rotation in stars that 
are significantly cooler than T e s = 6000 K. 

Considering stars with effec tive temperatu res in the range 
3400 K < r eff < 6000 K. IBarnes et al.1 J2005I) found a power- 
law dependence in their measurements of absolute shear AQ on 
stellar surface temperature. Their result is compared to the sam- 
ple of thi s work in Fig.|5] Absolute she ar follows the trend ex- 
pected by lKiiker & Riidigerl J2005h and lKitchatinov & RUd iger 
( 1999) with a stronger surface shear at higher effective tem- 
peratures. Although a large scatter is observed in the F-star 



measurements, they qualitatively follow this trend and connect 
to cooler stars at roughly the expected values. In addition, the 
rotation rate of the F-stars is indicated by symbol size, larger 
symbols displaying higher rotation rate. Stars exhibiting a shear 
in excess of the expected rate for their temperature tend to show 
very rapid rotation. Thus, the investigated F -stars agree qualita - 
tively with the temperature fit derived by Barn es et alJ (J2005). 
A more quantitative analysis is complicated by the large scatter 
among measure ments of absolute s hear, which is visible also 
in the sample of B arnes et alJ d2005l) . Again, it should be noted 
that the values plotted for the F-star sample only display stars 
for which signatures of differential rotation have been mea- 
sured. A large number of stars populate the region of weaker 
surface shear or rigid rotation (a = AQ = 0), and the tempera- 
ture law applies only to the strongest differential rotators. 

7. Conclusions 

The sample of stars with rotation laws measured from spectral 
broadening profiles with the FTM is constantly growing. In this 
work, 44 new observations of stars of spectral type F and later 
were added to the results from former publications. Currently, 
rotation laws have been analyzed in a homogeneous data set of 
147 stars of spectral type F and later, and in a second data set in 
78 stars of spectral type A. Among all these observations cov- 
ering the temperatur e range between 5600 K and 10000K (in- 
cluding A-stars from lReiners & RoverLl2004l) . 31 stars exhibit 
signatures of solar-like differential rotation. Only three of them 
are of spectral type A. In the HR-diagram, differential rotators 
appear near and on the cool side of the convection boundary. 
No differentially rotating star hotter than 7400 K is known, and 
it is obvious that the signatures of solar-like differential rota- 
tion are closely connected to the existence of deep convective 
envelopes. 

Most differential rotators can be found near to the ZAMS 
at young ages, but due to the limited sample and severe selec- 
tion effects, this needs confirmation from a less biased sam- 
ple. Particularly, a number of slower rotators with temperatures 
around T e ff = 6600 K are needed. 

Four differential rotators are very close to the convection 
boundary. All four show extraordinarily strong absolute shear 
and exhibit projected rotation velocities within 10% around 
v sin/ = 100km s . It is suggested that these stars form a 
group of objects in which rotation velocity and convection zone 
depth facilitate very strong absolute shear and that the mecha- 
nism causing the shear is different from the later F-type stars. 

Among the F-stars, differential rotation occurs in the whole 
range of temperatures and rotation rates. The sample of 147 
stars of spectral type F and later was investigated for the de- 
pendence of differential rotation on rotation and temperature. 
28 of them (19%) exhibit signatures of differential rotation. 
The distribution of differential rotation was approached with 
two different strategies: (i) investigation of the fraction of stars 
exhibiting differential rotation (a > 0); and (ii) analysis of the 
maximum a and maximum AQ as a function of rotation period 
and temperature. The first approach reflects the typical values 
of <a> and <AQ> at any given temperature and rotation rate, 
while the second focuses on the question how strong the abso- 
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lute shear can possibly be in such stars. Due to the large uncer- 
tainties in the measurements and the high minimum threshold 
in differential rotation, the mean values of a and AO are not 
particularly meaningful. Furthermore, it is not clear whether a 
smooth transition from stars exhibiting strong differential rota- 
tion to "rigidly" rotating stars (a < 0.05) or a distinction be- 
tween these two groups exists, and the mentioned approaches 
were preferred for the analysis. 

In the sample, hotter stars generally rotate more rapidly, and 
effects due to rotation velocity and temperature cannot be dis- 
entangled. The distribution of differential rotators depends on 
color and/or rotation rate, and the fraction of stars with differ- 
ential rotation (a > 0) increases with cooler temperature and/or 
slower rotation. It is not clear what this means to the mean 
shear, <AQ>. For example, it is not inconsistent with <AQ> 
being constant at all rotation rates. In this case, <a>—<A£l/Sl> 
would be smaller in more rapidly rotating stars, and thus a 
lower fraction of stars would exhibit differential rotation a 
above the observational threshold, as is observed. 

On the other hand, the maximum observed values of differ- 
ential rotation, a max , and of the absolute shear, A£2 max , do vary 
depending on the rotation rate. The strongest absolute shear of 
Afi w 1 radd~' is found at rotation periods between two and 
three days with significantly smaller values in slower rotators. 
The more rapidly rotating stars show a slight decrease in abso- 
lute shear as well, although the sparse data are also consistent 
with a plateau at AO * 0.7radd~' for 0.5 < P < 3d. A max- 
imum in differential rotatio n A£2 has recently been predicted 
bv lKiiker & Riidieerl ( Eool) for an F8-star, although of lesser 
strength and at a slower rotation rate, a difference that may in 
parts be due to their later spectral type. 

The investigated sample does not cover a wide range in ef- 
fective temperature since only very few late-type field stars ro- 
tate fast enough for the method applied. Although temperature 
is expected to strongly influence the strength of differential ro- 
tation, the large range in rotation rate and the connection be- 
tween rotation rate and temperature in the sample makes con- 
clusions about temperature effects insecure. The results were 
compared to differential rotation measurements in cooler stars 
and found in qualitative agreement with an extrapola t ion of 
the empirical temperature dependence that Bar nes et alj j2005) 
found when analyzing a sample of differentially rotation mea- 
surements done with DI. 

The implications of the relations discussed here for stel- 
lar magnetic activity and the nature of the dynamo working in 
F-type stars still remain unclear from an observational point 
of view. Naively, one would expect stronger magnetic activ- 
ity to occur in stars with stronger differential rotation among 
groups of comparable temperature or rotation rate. Those stars 
for which X-ray measurements are available, however, do not 
yet exhibit such a trend, but a meaningful investigation is ham- 
pered by the limited amount of data points available (especially 
for comparable temperature or rotation velocity). It has been 
shown that qualitative conclusions can be derived from the cur- 
rently available measurements of stellar rotation laws, but a 
more detailed investigation of the consequences on the dynamo 
operating in F-type stars has to wait until a statistically better- 
defined sample of stars is available. 
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Appendix A: The first direct comparison to 
Doppler Imaging 

Recently, IMarsden et alJ d2005h have presented Doppler 
Images of HD 307 938, a young active G dwarf in IC 2602. This 
star was also observed for this project in the FLAMES/UVES 
campaign, and I rep ort on the rotation law in Table lB.2l 
IMarsden et alJ d2005l) took a time series over four nights de- 
tecting spectroscopic variability. Their data is contaminated by 
a significant amount of sunlight reflected by the moon, which 
they have carefully removed before constructing Doppler 
Images. The result of Doppler Imaging is that HD 307 938 
has a cool polar cap extending down to ~ 60° latitude, and 
v sin/ = 92 + 0.5km , which excellently corresponds to 
v sin/ = 93.7 + 4.7kms~ I , as derived here (on comments about 
uncertainties in v sin i, see Sect.|6). They also report marginal 
differential rotation of AD. = 0.025 ± 0.015radd _1 with a l<x 
error. From this result, one expects a value of qilqx that is only 
marginally less than 1 .76, the value for rigid rotation and the 
best guess for the limb darkening parameter. The cool polar 
spot found on HD 307 938 influences qi/qi as well, as it en- 
larges qi/qi much more than the small deviation from rigid ro- 
tation does; q2/qi is thus expected to be larger than 1.76. The 
spectrum secured for the analysis presented here is not con- 
taminated by sunlight, and the broadening function derived is 
shown in Fig. IA.ll The profile does not reveal large asymmetry, 
although a spot may be visible around v ~ -30 km s _1 . The pro- 
fil e is fully consi s tent w ith the broadening functions presented 
in M arsden et al.1 (12005). who were able to show temporal vari- 
ations in the profile from lower quality data. The profile param- 
eter determined from FTM is qi/qi = 1.78 + 0.01, indicating 
no signs of solar-like differential rotation that is large enough 



to be detected with this method. However, the fact that q2/qi 
is slightly larger than 1.76 supports the idea that a cool spot 
occupies the polar caps. Thus, the finding of marginal differen- 
tial rotation and a cool polar spot is consistent with the result 
derived from FTM (besides the good consistency in v sin/). 
This is the first time that a direct comparison of the results is 
possible, since the majority of Doppler Imaging targets usually 
show spot signatures that are stronger than what can be dealt 
with using FTM. 

Appendix B: Tables of stars with measured 
differential rotation. 
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Table B.l. Field Stars with measured differential rotation. 



Star 


HR 


v sini 


6 v sin i 


qilqi 


Sqilqi 


AO, 


6A£l 




My 


PI sin i 






[kms^ 1 ] 


[kms- 1 ] 






rad d~' 


rad d~' 


K 


mag 


d 


HD 432 


21 


71.0 


3.6 


1.78 


0.03 


0.00 




6763 


0.96 


3.08 


HD 4089 


187 


23.5 


1.2 


1.82 


0.02 


0.00 




6161 


2.75 


5.01 


HD 4247 


197 


42.7 


2.1 


1.77 


0.04 


0.00 




6825 


3.08 


1.90 


HD 4757 


230 


93.8 


4.7 


1.73 


0.03 


0.00 




6706 


1.34 


2.00 


HD 6706 


329 


46.0 


2.3 


1.83 


0.06 


0.00 




6551 


2.77 


2.21 


HD 6903 


339 


87.7 


4.4 


1.69 


0.02 


0.54 


0.534 


6273 


3.38 


0.96 


HD 15524 


728 


59.8 


3.0 


1.81 


0.08 


0.00 




6592 


2.44 


1.95 


HD 17094 


813 


45.1 


2.3 


1.49 


0.02 


0.86 


0.231 


7141 


2.17 


2.50 


HD 17206 


818 


25.6 


1.3 


1.70 


0.02 


0.15 


0.184 


6371 


3.68 


2.78 


HD 18256 


869 


17.2 


0.9 


1.71 


0.04 


0.05 


0.121 


6332 


2.93 


5.91 


HD 22001 


1083 


13.2 


0.7 


1.83 


0.07 


0.00 




6601 


2.98 


6.88 


HD 22701 


1107 


55.0 


2.8 


1.83 


0.03 


0.00 




6610 


2.70 


1.88 


HD 23754 


1173 


13.8 


0.7 


1.87 


0.05 


0.00 




6518 


2.98 


6.75 


HD 24357 


1201 


65.8 


3.3 


1.82 


0.05 


0.00 




6895 


2.83 


1.35 


HD 25457 


1249 


18.0 


0.9 


1.71 


0.02 


0.09 


0.152 


6333 


4.07 


3.35 


HD 25621 


1257 


16.7 


0.8 


1.73 


0.03 


0.00 




6091 


2.26 


9.02 


HD 27459 


1356 


78.3 


3.9 


1.78 


0.05 


0.00 




7642 


2.15 


1.29 


HD 28677 


1432 


134.8 


6.7 


1.77 


0.03 


0.00 




6981 


2.76 


0.67 


HD 28704 


1434 


88.1 


4.4 


1.79 


0.03 


0.00 




6672 


2.43 


1.30 


HD 29875 


1502 


47.8 


2.4 


1.82 


0.06 


0.00 




7080 


3.39 


1.37 


HD 29992 


1503 


97.5 


4.9 


1.75 


0.04 


0.00 




6742 


2.64 


1.04 


HD 30034 


1507 


103.6 


5.2 


1.82 


0.06 


0.00 




7484 


2.34 


0.92 


HD 30652 


1543 


17.3 


0.9 


1.78 


0.03 


0.00 




6408 


3.59 


4.23 


HD 33167 


1668 


47.5 


2.4 


1.82 


0.03 


0.00 




6493 


2.10 


2.97 


HD 35296 


1780 


15.9 


0.8 


1.75 


0.02 


0.00 




6060 


4.11 


4.09 


HD 37147 


1905 


109.9 


5.5 


1.77 


0.07 


0.00 




7621 


2.25 


0.88 


HD 41074 


2132 


87.8 


4.4 


1.78 


0.06 


0.00 




6912 


2.45 


1.20 


HD 43386 


2241 


19.5 


1.0 


1.83 


0.03 


0.00 




6512 


3.50 


3.77 


HD 44497 


2287 


89.4 


4.5 


1.79 


0.03 


0.00 




7010 


1.79 


1.56 


HD 46273 


2384 


106.9 


5.3 


1.74 


0.05 


0.00 




6674 


2.10 


1.25 


HD 48737 


2484 


66.1 


3.3 


1.78 


0.04 


0.00 




6496 


2.32 


1.93 


HD 51199 


2590 


91.7 


4.6 


1.77 


0.04 


0.00 




6730 


2.06 


1.45 


HD 55052 


2706 


81.8 


4.1 


1.80 


0.04 


0.00 




6668 


0.63 


3.21 


HD 56986 


2777 


129.7 


6.5 


1.75 


0.08 


0.00 




6837 


1.95 


1.05 


HD 57927 


2816 


89.5 


4.5 


1.75 


0.02 


0.00 




6772 


1.40 


1.99 


HD 58579 


2837 


147.1 


7.4 


1.81 


0.03 


0.00 




7053 


1.69 


0.98 


HD 58946 


2852 


59.0 


3.0 


1.76 


0.07 


0.00 




6892 


2.79 


1.54 


HD 60111 


2887 


117.5 


5.9 


1.80 


0.05 


0.00 




7181 


3.01 


0.65 


HD 61035 


2926 


124.2 


6.2 


1.80 


0.03 


0.00 




6986 


3.09 


0.62 


HD 61110 


2930 


91.1 


4.6 


1.78 


0.04 


0.00 




6575 


1.51 


1.98 


HD 62952 


3015 


127.5 


6.4 


1.80 


0.08 


0.00 




6933 


1.53 


1.26 


HD 64685 


3087 


47.6 


2.4 


1.88 


0.03 


0.00 




6838 


3.04 


1.73 


HD 67483 


3184 


52.4 


4.1 


1.50 


0.07 


0.69 


0.279 


6209 


2.07 


3.01 


HD 69548 


3254 


53.9 


2.7 


1.77 


0.05 


0.00 




6705 


3.53 


1.27 


HD 70958 


3297 


45.5 


2.3 


1.74 


0.05 


0.00 




6230 


3.53 


1.76 


HD 70958 


3297 


46.1 


2.3 


1.75 


0.01 


0.00 




6230 


3.53 


1.74 


HD 72943 


3394 


56.8 


2.8 


1.70 


0.06 


0.17 


0.361 


6897 


1.81 


2.51 


HD 75486 


3505 


128.1 


6.4 


1.79 


0.02 


0.00 




6993 


1.22 


1.42 


HD 76143 


3537 


83.0 


4.2 


1.79 


0.03 


0.00 




6579 


1.58 


2.10 


HD 76582 


3565 


90.5 


4.5 


1.80 


0.05 


0.00 




7884 


2.58 


0.87 


HD 77370 


3598 


60.4 


3.0 


1.71 


0.03 


0.20 


0.408 


6609 


2.97 


1.51 


HD 77601 


3603 


140.7 


7.0 


1.82 


0.02 


0.00 




6421 


0.44 


2.21 


HD 79940 


3684 


117.2 


5.9 


1.78 


0.03 


0.00 




6397 


0.88 


2.18 


HD 8 1 997 


3759 


30.4 


1.5 


1.73 


0.01 


0.00 




6471 


3.28 


2.72 


HD 82554 


3795 


129.7 


6.5 


1.63 


0.05 


0.66 


0.460 


6272 


1.32 


1.68 


HD 83287 


3829 


102.5 


5.1 


1.77 


0.08 


0.00 




7815 


2.69 


0.74 


HD 83962 


3859 


140.3 


7.0 


1.72 


0.08 


0.00 




6507 


1.53 


1.30 


HD 84607 


3879 


93.1 


4.7 


1.79 


0.01 


0.00 




7000 


1.78 


1.51 


HD 88215 


3991 


97.5 


4.9 


1.78 


0.05 


0.00 




6823 


3.04 


0.85 


HD 89254 


4042 


63.5 


3.2 


1.81 


0.02 


0.00 




7173 


2.42 


1.57 
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Star 


rlK 


v sin i 


6 v sin i 


qilqx 


oqilqi 






T 

* eff 


My 


Pj sin / 






[Km s J 


|Km s J 






rad d _1 


rad d -1 




mag 


a 


HD 89449 


4054 


17.3 


1.7 


1.44 


0.04 


0.45 


0.127 


6398 


3.05 


5.44 


HD 89569 


4061 


12.2 


0.7 


1.57 


0.02 


0.21 


0.070 


6290 


3.10 


7.83 


HD 89571 


4062 


133.9 


6.7 


1.79 


0.07 


0.00 










HD 90089 


4084 


56.2 


2.8 


1.70 


0.09 


0.35 


0.987 


6674 


3.60 


1.19 


HD 90589 


4102 


51.6 


2.6 


1.93 


0.04 


0.00 




6794 


2.95 


1.68 


HD 96202 


4314 


93.4 


4.7 


1.77 


0.05 


0.00 




6747 


2.34 


1.25 


HD 99329 


4410 


137.9 


6.9 


1.73 


0.04 


0.00 




6990 


2.66 


0.68 


HD100563 


4455 


13.5 


0.7 


1.67 


0.04 


0.14 


0.134 


6489 


3.61 


5.24 


HD105452 


4623 


23.5 


1.2 


1.59 


0.02 


0.40 


0.151 


6839 


2.98 


3.60 


HD 106022 


4642 


77.2 


3.9 


1.86 


0.06 


0.00 




6651 


2.28 


1.60 


HD 107326 


4694 


132.2 


6.6 


1.80 


0.05 


0.00 




7107 


2.08 


0.90 


HD 108722 


4753 


97.0 


4.8 


1.78 


0.04 


0.00 




6490 


1.71 


1.74 


HD109085 


4775 


60.0 


3.0 


1.75 


0.05 


0.00 




6813 


3.13 


1.32 


HD109141 


4776 


135.7 


6.8 


1.79 


0.03 


0.00 




6881 


2.84 


0.66 


HD1 10385 


4827 


105.2 


5.3 


1.75 


0.01 


0.00 




6717 


1.71 


1.49 


HD1 10834 


4843 


133.3 


6.7 


1.73 


0.06 


0.00 




6244 


0.99 


1.93 


HD111812 


4883 


63.0 


3.2 


1.78 


0.05 


0.00 




5623 


2.87 


2.22 


HD 11 2429 


4916 


119.6 


6.0 


1.77 


0.03 


0.00 




7126 


2.97 


0.66 


HD 114378 


4968 


19.9 


1.0 


1.76 


0.01 


0.00 




6324 


3.82 


3.40 


HD115810 


5025 


99.2 


5.0 


1.82 


0.06 


0.00 




7185 


1.95 


1.24 


HD 116568 


5050 


36.8 


1.8 


1.73 


0.01 


0.00 




6485 


3.20 


2.32 


HD118889 


5138 


140.6 


7.0 


1.82 


0.06 


0.00 




6951 


2.40 


0.76 


HD1 19756 


5168 


63.9 


3.2 


1.78 


0.03 


0.00 




6809 


3.06 


1.29 


HD120136 


5185 


15.6 


1.0 


1.57 


0.04 


0.31 


0.134 


6437 


3.38 


5.13 


HD121370 


5235 


13.5 


1.3 


1.46 


0.03 


0.21 


0.059 


6024 


2.36 


0.90 


HD 122066 


5257 


40.6 


2.0 


1.81 


0.01 


0.00 




6395 


2.18 


3.45 


HD124780 


5337 


70.7 


3.5 


1.80 


0.04 


0.00 




7204 


2.21 


1.54 


HD124850 


5338 


15.0 


0.8 


1.91 


0.04 


0.00 




6075 


2.85 


7.71 


HD 127739 


5434 


55.8 


2.8 


1.81 


0.01 


0.00 




6787 


2.20 


2.20 


HD127821 


5436 


55.6 


2.8 


1.75 


0.04 


0.00 




6601 


3.76 


1.14 


HD129153 


5473 


105.7 


5.3 


1.78 


0.06 


0.00 




7693 


2.70 


0.73 


HD 129502 


5487 


47.0 


2.4 


1.80 


0.03 


0.00 




6695 


2.94 


1.91 


HD 129926 


5497 


112.5 


5.6 


1.74 


0.03 


0.00 




6048 


3.91 


0.64 


HD132052 


5570 


113.2 


5.7 


1.77 


0.06 


0.00 




6964 


2.19 


1.03 


HD136359 


5700 


20.3 


1.0 


1.78 


0.01 


0.00 




6296 


3.02 


4.86 


HD136751 


5716 


72.7 


3.6 


1.77 


0.02 


0.00 




6810 


2.31 


1.60 


HD138917 


5788 


85.8 


5.4 


1.99 


0.06 


0.00 










HD 139225 


5804 


104.4 


5.2 


1.76 


0.03 


0.00 




6960 


2.50 


0.97 


HD 139664 


5825 


71.6 


3.6 


1.77 


0.05 


0.00 




6681 


3.57 


0.94 


HD 142908 


5936 


75.8 


3.8 


1.80 


0.05 


0.00 




6849 


2.33 


1.50 


HD 143466 


5960 


141.3 


7.1 


1.77 


0.05 


0.00 




7235 


2.20 


0.77 


HD 147365 


6091 


72.5 


3.6 


1.77 


0.06 


0.00 




6657 


3.46 


0.99 


HD 147449 


6093 


76.4 


3.8 


1.79 


0.04 


0.00 




6973 


2.70 


1.21 


HD 147449 


6093 


77.1 


3.9 


1.80 


0.03 


0.00 




6973 


2.70 


1.20 


HD 148048 


6116 


84.8 


4.2 


1.78 


0.05 


0.00 




6731 


2.43 


1.32 


HD150557 


6205 


61.8 


3.1 


1.83 


0.04 


0.00 




6959 


2.01 


2.06 


HD151613 


6237 


47.5 


2.4 


1.84 


0.07 


0.00 




6630 


2.80 


2.06 


HD155103 


6377 


57.9 


2.9 


1.81 


0.05 


0.00 




7150 


2.55 


1.63 


HD156295 


6421 


107.4 


5.4 


1.79 


0.07 


0.00 




7818 


2.56 


0.75 


HD160915 


6595 


12.4 


0.6 


1.60 


0.03 


0.23 


0.101 


6356 


3.58 


6.07 


HD164259 


6710 


69.3 


3.5 


1.75 


0.04 


0.00 




6704 


2.70 


1.44 


HD165373 


6754 


79.9 


4.0 


1.74 


0.03 


0.00 




6976 


2.08 


1.54 


HD171834 


6987 


71.3 


3.6 


1.77 


0.04 


0.00 




6622 


2.50 


1.58 


HD173417 


7044 


53.9 


2.7 


1.83 


0.04 


0.00 




6780 


1.81 


2.74 


HD173667 


7061 


18.0 


2.0 


1.40 


0.02 


0.44 


0.104 


6363 


2.78 


5.99 


HD175317 


7126 


17.1 


0.9 


1.58 


0.04 


0.30 


0.136 


6563 


3.11 


5.07 


HD 175 824 


7154 


53.7 


2.7 


1.69 


0.03 


0.15 


0.179 


6232 


1.72 


3.43 


HD 182640 


7377 


87.3 


4.4 


1.67 


0.06 


0.62 


0.753 


7016 


2.46 


1.17 


HD185124 


7460 


87.0 


4.4 


1.71 


0.03 


0.30 


0.604 


6680 


2.98 


1.02 


HD 186005 


7489 


149.9 


7.5 


1.79 


0.04 


0.00 




6988 


1.68 


0.98 
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Star 


rlK 


v sin i 


S v sin i 


qilqx 


oqilqi 






1 eff 


My 


Pj sin / 






[Km s J 


n™ „ — 1 1 

[Km s J 






rad d _1 


,1 J-l 

rad a 


1/ 

K 


mag 


a 


HD187532 


7553 


77.5 


3.9 


1.75 


0.03 


0.00 




6788 


3.34 


0.94 


HD 189245 


7631 


72.6 


3.6 


1.74 


0.03 


0.00 




6259 


4.06 


0.85 


HD 190004 


7657 


136.1 


6.8 


1.80 


0.04 


0.00 




6974 


2.47 


0.76 


HD 197692 


7936 


41.7 


2.1 


1.62 


0.02 


0.65 


0.286 


6587 


3.33 


1.86 


HD 199260 


8013 


13.7 


0.7 


1.79 


0.03 


0.00 




6213 


4.18 


4.36 


HD201636 


8099 


58.8 


2.9 


1.86 


0.03 


0.00 




6700 


2.20 


2.15 


HD203925 


8198 


70.7 


3.5 


1.80 


0.03 


0.00 




6845 


1.15 


2.77 


HD205289 


8245 


57.5 


2.9 


1.74 


0.03 


0.00 




6525 


3.24 


1.44 


HD206043 


8276 


134.0 


6.7 


1.74 


0.04 


0.00 




7092 


2.87 


0.62 


HD207958 


8351 


69.3 


3.5 


1.84 


0.03 


0.00 




6747 


2.97 


1.26 


HD2 10302 


8447 


13.6 


0.7 


1.72 


0.04 


0.00 




6465 


3.52 


5.46 


HD210459 


8454 


147.4 


7.4 


1.76 


0.03 


0.00 




6376 


0.21 


2.89 


HD2 13051 


8558 


48.4 


2.4 


1.85 


0.08 


0.00 










HD213845 


8592 


35.7 


1.8 


1.71 


0.02 


0.14 


0.234 


6551 


3.35 


2.18 


HD2 19693 


8859 


19.9 


1.0 


1.81 


0.01 


0.00 




6461 


2.82 


5.14 


HD220657 


8905 


73.4 


3.7 


1.74 


0.03 


0.00 




5801 


2.41 


2.15 



Table B.2. Stars in cluster fields, FLAMES/UVES observations 



Star 


Cluster 


v sin i 


6 v sin i 


qilqi 


Sqilq\ 


AO, 


SAQ, 




M v 


PI sin 






[kms- 1 ] 


[kms" 1 ] 






rad d- 1 


radd- 1 


K 


mag 


d 


NGC 6475 69 


NGC 6475 


94.9 


4.7 


1.76 


0.01 


0.00 




6822 


3.27 


0.75 


NGC 6475 41 


NGC 6475 


80.1 


4.0 


1.71 


0.01 


0.41 


0.497 


6346 


3.86 


0.76 


BD+20 2161 


Praesepe 


72.7 


3.6 


1.76 


0.01 


0.00 




6864 


3.23 


0.99 


BD+20 2170 


Praesepe 


94.4 


4.7 


1.71 


0.01 


0.46 


0.557 


6484 


3.69 


0.68 


CI* NGC 2632 KW 230 


Praesepe 


78.9 


3.9 


1.75 


0.01 


0.00 




6628 


3.51 


0.85 


CI* IC 2391 L 33 


IC2391 


81.1 


4.1 


1.74 


0.01 


0.00 




6414 


3.78 


0.77 


HD 307938 


IC 2602 


93.7 


4.7 


1.78 


0.01 


0.00 




5743 


4.64 


0.51 


HD 308012 


IC2602 


45.7 


2.3 


1.70 


0.01 


0.36 


0.306 


6147 


4.12 


1.23 


CI* IC 4665 V 69 


IC4665 


45.1 


2.3 


1.74 


0.01 


0.00 




6180 


4.08 


1.26 


CI* IC 4665 V 102 


IC4665 


105.0 


12.0 


1.39 


0.01 


3.62 


0.776 


7136 


2.93 


0.73 


CI* IC 4665 V 97 


IC4665 


52.0 


2.6 


1.90 


0.01 


0.00 




6628 


3.51 


1.29 



